Multi-patch model for transport properties of cuprate superconductors 
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A number of normal state transport properties of cuprate superconductors are analyzed in detail 
using the Boltzmann equation. The momentum dependence of the electronic structure and the 
strong momentum anisotropy of the electronic scattering are included in a phenomenological way 
via a multi-patch model. The Brillouin zone and the Fermi surface are divided in regions where 
scattering between the electrons is strong and the Fermi velocity is low (hot patches) and in regions 
where the scattering is weak and the Fermi velocity is large (cold patches). We present several 
motivations for this phenomenology starting from various microscopic approaches. A solution of 
the Boltzmann equation in the case of N patches is obtained and an expression for the distribution 
function away from equilibrium is given. Within this framework, and limiting our analysis to 
the two patches case, the temperature dependence of resistivity, thermoelectric power. Hall angle, 
magnetoresistance and thermal Hall conductivity are studied in a systematic way analyzing the 
role of the patch geometry and the temperature dependence of the scattering rates. In the case 
of Bi-based cuprates, using ARPES data for the electronic structure, and assuming an inter-patch 
scattering between hot and cold states with a linear temperature dependence, a reasonable agreement 
with the available experiments is obtained. 
PACS numbers: 72.10.Di, 74.25. Fy, 71.10.Ay 
I. INTRODUCTION 



The normal state transport properties of cuprate su- 
perconductors have attracted enormous attention. In the 
underdoped regime a pseudogap appears in the excitation 
spectrum of the metallic state above the superconduct- 
ing critical temperature Tc and below a doping dependent 
crossover temperature T*. At optimum doping the T* al- 
most coincides with Tc and the pseudogap region of the 
phase diagram, if present, is very narrow. On the other 
hand the metallic state of cuprates at optimum doping, 
away from the pseudogap region, still strongly deviates 
from what is observed in simple metals. At optimum 
doping the in-plane DC-resistivity is liiWgf in tempera- 
ture from Tc to very liigh temperatureaj'cl, the thermo- 
electric power is lineam, the cotangent of the Hall angle 



displays a T^- 
based cupratei 



ependence, (7 ~ 2 in La-based and Y- 



'0, with deviations at low temperature, of 
the order of 2Tc, 1.60 < 7 < 2 in Bi-hased cupratesQ'tO, 
where increasing the doping decreases the exponent 7), 
the magnetoresistagce has approximately a T~" depen- 
dence, with a ~ 40 and the thermal Hall conductivi 
has approximately a T~^ dependence, with /3 ~ 1.21 
Increasing the doping, in the overdoped regime, the con- 
ventional Fermi liquid character of the metallic state is 
almost recovered. 

Theoretical approaches to the problem of transport 
properties of cuprates can be classified in: (z) ap- 
proaches where in-plane transport properties are ana- 
lyzed in terms of two scattering times, one associated 
to the response to an electric field (determining DC- 
conductivity) and the other one associated to the re- 
sponse to a magnetic field, (determining Hall conductivity 
and magnetoresist|aiice)ll3; (ii) approaches based on the 
Boltzmann theoryllj, where the scattering time is mo- 
mentum dependent and the Fermi surface can be divided 



in hot regions (around the M points of the Brillouin zone 
(BZ)) corresponding to strong scattering between quasi- 
particles (short scattering time) and low Fermi velocity 
and cold regions (around the nodal points located along 
the diagonals FY and FX of the BZ), corresponding to 
weak scattering (large scatteri|M ptime) and large Fermi 
velocity. In the various modelal^il3 proposed within the 
class (a), the temperature dependence of the scattering 
time in the hot regions strongly deviates from the 
behavior of simple metals, while recovering this conven- 
tional behavior in the cold regions. Hot/cold regions (or 
spots) models are able to capture some anomalous prop- 
erties of cuprates, but a general consensus and a sys- 
tematic analysis of the full set of electric and thermal 
transport properties is lacking in the literature. 

In this paper we study in a systematic way the nor- 
mal state transport properties of cuprate superconduc- 
tors within the second approach. We introduce a new 
parameterization of the scattering matrix in the Boltz- 
mann equation (BE) for the quasiparticle distribution 
function via a multi-patch model, motivated by the strong 
momentum dependence of the electronic properties ob- 
served in the cuprates in angle resolved photoemission 
spectroscopy (ARPES). The BE is solved in the case of 
N patches and an expression for the perturbed quasi- 
particle distribution function is obtained. Our analysis 
is then limited for simplicity to the case of two patches, 
leaving the multi-patch problem for future investigations. 
Within the two-patch model the BZ and the Fermi sur- 
face (FS) are divided into two regions, the hot regions cor- 
responding to hot patches on the FS and Ihe cold regions 
corresponding to cold patches on the FSt£l. Assuming a 
temperature dependence for the scattering amplitude 
in the cold region, a T dependence for the inter-patch 
(hot-cold) scattering and a constant temperature depen- 
dence in the hot region, a reasonable description of the 
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available experimental data for _Bi-based cuprates is ob- 
tained. 

ARPES experiments, mainly performed in Bi22l2 
compounds, give a qualitative justification to the divi- 
sion of the BZ and of the FS in cold and hot regions. 
ARPES clearly shows that the line-shape (Atiip spectral 
function is strongly momentum dependentt3'ta. Around 
the M points of the Brillouin zone (±7r,0); (0,±7r) the 
spectral function is very broad (thsJine-width is of the 
order of 0.2 - O.SeV at T = lOOidla) and a quasiparti- 
cle peak cannot be easily distinguished; the states around 
the M points are therefore almost incoherent (as localized 
states) and a very strong scattering mechanism is at the 
origin of the broad line-shape. Moreover the electronic 
band dispersion has saddle points located at the M points 
and this originates the van Hove singularity in the density 
of states. The band dispersion along the AIY direction is 
very narrow (Ae « bOmeV foe. optimally doped _Bz2212) 
and the Fermi velocity is lowllj. These states correspond 
therefore to hot states. The incoherent behavior and the 
associated line-shape is also temperature independent in 
the wide range of temperature between and 300K. On 
the other hand, around the nodal points of the BZ located 
along the rY{X) diagonals, the spectral function has a 
well pronounced quasiparticle peak (thSj line- width is of 
the order of 0.05 - O.leV at T = lOOidla) and the wave- 
vector dispersion of this peak together with the tem- 
perature dependence of the peak width can be followed. 
Valla et al. found a linear temperature and frequency de- 
pendenee of the peak width for states around the nodal 
pointsE2l. Moreover the band dispersion along TY{X) is 
wide (Ae « AOQmeV for optimally doped i3j2212) and 
the Fermi velocity is high. The ratio of the Femii ve- 
locities in the two regions is VF{rY)/vF{M) ~ 3cil. The 
quasiparticle states around the nodal points are therefore 
coherent (delocalized states) and a scattering mechanism 
with weaker intensity (even with unconventional nature) 
is at the origin of the line-shape behavior. These states 
correspond to cold states. 

The physics of cuprates is very rich. Starting from the 
low doping insulating phase, the holes added to the Cu02 
planes of cuprates through out-of-plane doping will be- 
come metallic and segments of a FS appear. The holes 
move in an antiferromagnetic background and two holes 
with opposite spins on the same lattice site experience 
a strong Hubbard repulsion. The conducting holes are 
therefore strongly correlated and a possible description of 
the electronic properties of the metallic state of cuprates 
has been formulated in terms of the two dimensional 
Hubbard model, eventually with the inclusion of other 
terms in the Hamiltonian (extended Hubbard model) to 
take into account short range attraction due to phonons 
and/or long range Coulomb interaction which can be rel- 
evant at high doping. First insights in the phase diagram 
(temperature vs. doping) of the (extended) Hubbard 
model indicate the presence of several electronic instabil- 
ities arising from the competition between the different 
degrees of freedom present in the Hamiltonian, such as 



the kinetic term (delocalizing term) and the short range 
terms (localizing terms). The most relevant phases, to- 
ward which the system is unstable, are the antiferromag- 
netic insulating phase, phase separation in macroscopic 
regions with low and high density of holes, spin and 
charge ordering in the form of stripes, and finally super- 
conductivity. The metallic state of cuprates, in particular 
in the underdoped region of the phase diagram, can be 
close to one of these electronic instabilities. Therefore, 
the properties of the metallic state of cuprates can be 
strongly connected to the presence of several competing 
interactions which arise nearby the electronic instabilities 
mentioned above. 

There are two possible microscopic origins for the mo- 
mentum differentiation of the BZ, one connected to elec- 
tronic scattering mediated by spin, charge or pair fluctu- 
ations and another associated with proximity to a Mott 
transition. 

The first possibility for the momentum differentia- 
tion of the BZ is based on its connection to electronic 
scattering due to spin fluctuation^. r-Superconducting 
fluctuationsllil and charge instabilitiescS have also been 
involved. These fluctuations can mediate the electron- 
electron interaction in both the particle-hole and particle- 
particle channels and hence they can determine strong 
deviations in the properties of the metallic state, such 
as the line-shape of the ARPES spectra and the shape 
of the FS respect to a conventional Fermi liquid picture. 
The various fluctuations have a speciflc momentum and 
frequency dependence and the propagators are peaked at 
different critical momenta qc, depending from the under- 
lying instability: for phase separation qc t^oEj, for an- 
tiferromagnetic insulatprjOc = Q = (tt, 7r)E3, for charge 
ordering qc = qstripeEju, where Qstripe = ^n/Xstripe 
[Xstripe is the pcriodicity of the charge modulation). 
Therefore, the critical fluctuations couple electrons with 
momenta only inside particular segments of FS, and the 
specific geometry of the strongly coupled (and hence of 
the weakly coupled) segments is determined by the inter- 
play between the shape of the FS and qc. This leads to 
the phenomenology of the two- or multi-patch model for 
the electronic scattering in the metallic state of cuprates. 

The second scenario for momentum differentiation in- 
vokes the proximity to a Mott transition. It builds on our 
recent understanding of this pheaiomena within dynam- 
ical mean field theory|-(pMFT)o and its extension: the 
two impurity methodEJ'Ea, the Bethe-PeieclS|-dustcrE3, 
the dynamical cluster approximation (DCA)E3"cii«id| jtlie 
cellular dynamical mean field theory (C-DMFT)E3t!3^l!a. 

DMFT allows a microscopic description of the strongly 
correlated state near the Mott transition. There is a 
temperature scale reminiscent of the Kondo temperature 
Tfj-, such that for T ^ the quasiparticles are Fermi 
liquid like, while for T > Tk the single particle exci- 
tations become incoherent and the transport properties 
are non Fermi liquid like. In single site DMFT studies, 
the two regimes were obtained by varying the tempera- 
ture and the strength of the local Hubbard repulsion [/, 
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and by construction they occur uniformly in the BZ. Re- 
laxing the constraint of momentum independent selfen- 
ergy by using the C-DMFT, one envisions that different 
patches of BZ have different Kondo temperatures, lead- 
ing naturally to the multi-patch model presented here. 
According to this view, therefore, the microscopic origin 
of the multi-patch division of the BZ and of the FS of 
cuprate superconductors can be associated to the pres- 
ence of a nearby transition from a (non conventional) 
metallic state to a Mott antiferromagnetic insulator. A 
detailed microscopic derivation of the many-patch model 
is in progress, using the C-DMFT. 

These views on momentum space differentiation may 
be complementary, and there are already strong hints 
from numerical calculations by Onoda and ImadaE3, that 
they occur in the Hubbard model. 

The plan of the paper is the following. In Section || 
we recall the BE in the linearized form and we introduce 
the multi-patch model for the the scattering operator. 
The scattering operator of the BE is projected on the 
patches and temperature dependences of its coefficients 
are assigned. A set of smooth functions is introduced to 
permit a continuous transition between hot and cold re- 
gions. A solution of the BE in the case of N patches is 
given in terms of the perturbed quasiparticle distribution 
function. In Section II] the results for the normal state 



transport properties obtained by the two-patch model are 
presented in a systematic way. The temperature depen- 
dences of resistivity, thermoelectric power, cotangent of 
the Hall angle, magnetoresistance and thermal Hall con- 
ductivity are reported. The various transport properties 
are studied for different set of parameters and the rele- 
vant hot /cold patch is associated to every quantity. Our 
model is applied to the normal state transport properties 
of _Bz-based cuprates {Bi2212 and Bi2201) and, taking 
ARPES data as input for the electronic structure, we find 
a reasonable agreement with the available experimental 
data. Discussions and conclusions are given in Section 
IV. 



II. THE BOLTZMANN EQUATION AND THE 
MULTI-PATCH MODEL 

ARPES experiments show that optimally and over- 
doped cuprates have a large FS in the normal state and 
well defined quasiparticles exist in a sizeable (cold) re- 
gion of the BZ around the rY{X) directions. The ex- 
istence of a FS and quasiparticles makes the treatment 
within the framework of the BE possible. We start our 
analysis introducing the linearized BE. As we are inter- 
ested in electrical, heat and Hall transport properties, 
we consider three terms in the BE, the terms including 
the electric field E and the temperature gradient WrT/T 
(driving terms) and the term including the magnetic field 
B (bending term). We consider here uniform electric and 
magnetic field. In this case the linearized BE has the fol- 
lowing form: 



dt 



9k 



9ek 



EkVk 



eEvk 



VrT 
T 

ft^t"'^"^]- air 



= Ck (1) 



where 5k = ffk + {-df^)/{dew)J2k'h,k'gk' takes the 
interaction between quasiparticles /k,k' into account 
(Fermi liquid corrections)^] and 5k is the departure from 
the equilibrium distribution function Vk is the group 
velocity of the quasiparticles. The scattering operator 
Ck has the form Ck = [Ck,k'5k' - Ck^k'9k], where 
C'k.k' is the scattering matrix, describing the scatter- 
ing of quasiparticles on an effective bosonic mode or 
impurity centers with the first term describing scatter- 
ing "in" to the state k and the second term describ- 
ing scattering "out" of the state k. The relaxation 
time Tk for the state k is defined as l/rk = X^k' C'k,k'- 
Within the conventional microscopic approach to the 
transport properties of a Fermi liquid, the transport 
equation for the distribution function can be written 
in terms of the four-point vertex part r(k, k', q, lj)EJ. 
The scattering matrix Ck,k' is then identified with the 
q = 0, a; = limit of the irreducible part of the vertex 
Ck,k' =-i(zV2)r(^Hk,k',q = 0,u; = 0), where is the 
mass renormalization and F*-^-* is the irreducible vertex. 
Below we will use this identification to give a first micro- 
scopic justification to our choice of Ck,k'- In the steady- 
state case we can replace 5k — > 5k as every term in Eq. 
(|l|) is expressed in terms of g^- Therefore the knowledge 
of the form of /k,k' is not required in the steady-state 
case. Frequency-dependent transport processes require 
both 5k and 5k, and hence the knowledge of /k,k' be- 
comes important. In the steady-state case the number of 
quasiparticles gk contributing to transport is given by 

" e 1 1 

— (vk X B) • Vk + — 5k - Ck,k'5k' = (2) 



evk • E - EkVk • 



dek 



Defining the operator ^k,k' as 



Ak.k' = 



— + ^(vk X B) 

Tk he 



3k,k' 



Ck.l 



(3) 



allows us to write the l.h.s. of Eq. ^ in the form 

J2k' ^k,k'5k' and hence the inverse of Ak,k' is required 
to solve Eq. (||). Considering weak magnetic fields, 
the bending term containing the magnetic field in the 
BE can be treated as a small perturbation of the trans- 
port process. We split ^k,k' in two parts, Ak,k' = 
Kk.k' + ^^kk'' with a magnetic field independent part 
^k,k' = (l/7'k)<5k.k' — C'k.k' and a part that contains 
the magnetic field M^^, = [(e//ic)(vk x B) • Vk](5k,k'- 

A perturbative expansion of Ak,k' in powers of B allows 
us to write the inverse of this operator as 
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A 



k,k' 



if, 



k,k' 



^k,L^k^. 



-^k.ki^'^^kSkJ-f^kj, kk^"kk,kl-f^kl k' 



k;. 

Ml 



kJ^'kJ.k 

B 



(4) 



K- 



with a summation over repeated indexes. Depending on 
the quantity of interest, we get contribution from the 
different terms in this expansion. The first term gives 
the leading order contribution to the DC and thermal 
conductivity, the second term to the (thermal) Hall- 
conductivity and the third term to the magnetoresis- 
tance. It follows from Eq. (||) that the number of parti- 
cles (7k contributing to transport is given by 



5k 



1 

k,k' 



evk' • E - ek'Vk' 



T 



9ek' 



(5) 



This equation is the starting point of our analysis. Elec- 
trical and thermal conductivities are derived using gk 
obtained from Eq. (||). 

The transport properties are separable in electric and 
thermal properties. All possible currents are given by 



Je 
JQ 



a_ S 
S R 



E 

VT 
T 



(6) 



where je is the electric current, jq is the thermal current, 
a is the electrical conductivity tensor, R is the thermal 
conductivity tensor and S is the thermopower tensor. 
Note that Eq. (^) contains a symmetric matrix using 
VT/T as driving term for thermal gradients. The cur- 
rents defined in Eq. (j^) can also be expressed in terms of 
quasiparticles. The electric and the thermal currents are 
given by je = e^^k^k^k and jg = X^k^kCk^k- (When 
frequency dependence is considered, ^k has to be replaced 
with 5k in the expression of the currents given above.) 
The tensors defined in Eq. (0) are given by: 



III' 



2e^E<^k,k'<' 

k,k' 



k,k' 



k.k'^k'Wk' 



-2e5^<ekik,l' 

k,k' 



9ek' 



9ek' 



9ek' 



(7) 



(8) 



(9) 



with the inverse of the operator Ak,k' given in Eq. (jj). 
The factor 2 in the expressions above takes the spin de- 
generacy into account. 

The last unknown quantity in Eq. (|^) is the scatter- 
ing matrix Ck,k'- The scattering matrix is connected 
via a frequency integration to the spectral function of 
the effective bosonic mode exchanged in the electronic 
scattering. As discussed above, ARPES suggests that 
the bosonic mode is strongly momentum dependent and 



divides the BZ in hot and cold regions. In order to in- 
clude in the scattering matrix a non trivial momentum 
dependence and the cold/hot division of the BZ, while 
maintaining a simple solution of the BE, a possibility 
is to expand Ck,k' with respect to a basis of functions 
which are able to select the various regions of the BZ, 
according to the momentum dependence of the effective 
bosonic mode. Therefore, the scattering matrix Ck,k' can 
be written as 



N 



a,k' = J2 «..*.(k)<i>,(k'), 



(10) 



where ^^(k) is a function which is equal to one inside 
the i-th patch of the BZ and zero outside, and it inter- 
polates continuously between these two values; aij is the 
amplitude of the scattering between the i-th and the j-th 
region of the BZ, it is in general temperature dependent 
and {{aij)) is a symmetric matrix with real elements; N 
is the total number of regions of the BZ required to in- 
clude properly the main effects of the anisotropy of the 
scattering and its interplay with the shape of the FS. 

One possible microscopic motivation for the form 
of the scattering matrix given in Eq. ( |l0| ) is 
based on the connection between Ck.k' and the four- 
point vertex F discussed above Eq. (|^) and on 
the C-DMFT. The four-point vertex can in princi- 
ple be evaluated by C-DMFT and the general struc- 
ture of the irreducible vertex will be F^^^ (k, k', q, a;) = 
(a;)</.l(k)00(k')</>t(k + q)0,(k' + q). 



V F'^ 
Once the patch function is defined as 



= (",/3) 



(k) 



(/)Jj(k)0^(k), Eq. (10) for Ck,k' can be recovered. 
The coefficient Uij entering Eq. ( |lO| ) are then de- 
fined in terms of the four-point vertex of the cluster as 
a^=|^a,0),J=i',,S) =i{z^/2)T''^f^^^g{uj = 0). 

The form of the scattering matrix given in Eq. ( p^ ) 
permits an analytical solution of the linearized BE for 
gk defined above. Considering the electric field and the 
thermal gradient as external perturbations, Eq. (|^) gives 



5k = Tk 



evkE - ekVk 



VrT 



T 



N 



U)+rk^a.,a>.(k). (11) 



■E 



eVkE - EkVk 



VrT 

T 



N 



which reduces the solution of the BE to a matrix inversion 
problem for the elements of a matrix Tji given by 



N 



Tjl = Sjl - E ""J E *m(k)Tk$/(k). 



(12) 



m— 1 



From Eq. (|Tl|) and Eq. (|T|) we can obtain the inverse 
of the matrix i^k.k' , defined above Eq. (^) , as follow 



4 



N N 

1=1 

(13) 

Once wc know the explicit expression for the operator 
k'' ^^'^ ^t)le to evaluate the expansion of the op- 
erator A^]^, for weak magnetic fields given in Eq. (|^). 

Inserting the operator A^^\^, in Eq. (^), we can evaluate 
the solution of the linearized BE, given by (;k, in the pres- 
ence of an electric field, a thermal gradient and a weak 
magnetic field, gk can be then used to evaluate all the 
currents given above Eq. (Q). All the response functions 
given in Eqs. ||, ^ are obtained inserting directly the 
expression for . 

In the following we consider a minimal realization of 
this multi-patch approximation limiting our analysis to 
N=2 patches. The two-patch model permits the distinc- 
tion between cold and hot regions on the BZ and the 
FS and it is suitable to describe the main properties of 
the scattering between electrons and an effective mode 
peaked at large momenta {e.g. an antifcrromagnetic spin 
fluctuation peaked at Q = (tt; tt)). On the other hand, to 
include the effect of the (small momenta) forward scatter- 
ing, important for transport mainly in the cold region, a 
larger number of patches is required (at least N=5), and 
we deserve this case for further investigation. 

In the case N=2, the two-patch division of the BZ and 
of the Fermi surface is realized by introducing two func- 
tions <I>k and ^I^k: <&k describes the cold region and ^fk 
the hot region as indicated in Fig. |^. 




FIG. 1. The first quadrant of the BZ is splitted into cold 
and hot regions. Cold regions are described by a function <l?k 
and hot regions by a function vl'k- The angle 6 parameterizes 
the size of the hot regions. 

The cold region is described by four boundaries as 
shown in Fig. [l[ The angle 9, defined with respect 
to the Y — (tTjTt) point of the BZ, parameterizes the 
size of the hot regions (and hence of the cold regions). 



The reason to split the whole BZ (and not only the FS) 
into a cold and a hot part is that the term df^/de in 
Eq. (||) doesn't restrict the sum over k on the FS any 
more when we increase size-ably the temperature, study- 
ing T-dependent properties. As the leading contribu- 
tion to the magnetoresistance is given by the third term 
in Eq. (^), which contains two derivatives with respect 
to kxjky (given by repeated two times), the mag- 
netoresistance diverges in the case of a discontinuous 
change between the two patches. Therefore we intro- 
duce functions ^k and 4'k varying smoothly between the 
two regions. We use hyperbolic tangents to describe the 
smooth change between the two regions and a parameter 
w is introduced to describe the width of the transition 
region. In the limit w — > a step function is recov- 
ered, d{x) — limt„^o(l + ta,nh{x/w))/2. Finally, we de- 
fine two slopes, mi and m2, and two offsets, ti and t2, 
given by mi = cot(0) > 0, m2 — tan(6') > 0, mi > m2, 
ti = 7r(l - tan(6')) > and t2 = n{l ~ cot(6')) < 0. In 
the case of smooth functions, the cold region in the first 
quadrant of the BZ (see Fig. n^) is described by: 



$k=n**w 



(14) 



with 



*i(k) 
$2(k) 
$3(k) 
*4(k) 



1 — tanh 



tanh 



w 



m2kx 



w 



1 — tanh 



tanh 



ky — m2kx — ti 



mikx - <2 



The smooth change between the two patches is shown in 
Fig. I 



1*. 



«I>. 



COLD 




FIG. 2. Smooth change between the two patches. The pa- 
rameter w describes the width of the transition region between 
the two patches and it allows to study the effects of this region 
on transport properties. 

The hot region is described by ^'k = 1 — 'I'k- Note that 
^^k vanishes only asymptotically in the cold region and 
vice versa. 

Following the phenomenological approach described 
above, the scattering matrix can be written consisting 
of different scattering mechanisms in different patches. 
In the case of the two-patch model, we introduce three 
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parameters that describe the possible scattering, scatter- 
ing inside the cold or inside the hot region and scattering 
between the two (hot/cold) regions, 

Ck,k' = a$k*k' + ^'^k^'k' + c[$k*k' + *k*k']. (15) 

The scattering matrix Ck,k' given in Eq. (|l5| ) is written 
as a sum of terms which have a separate dependence from 
k and k'. This approximation permits an analytic solu- 
tion of the BE for ^k, while having a non trivial momen- 
tum dependence and symmetry properties of the scatter- 
ing process. With the symmetric scattering matrix given 
in Eq. (p^), we obtain a momentum dependent relax- 
ation time Tk, 1/Tk = J2k' C'k.k', 



1 



Tk 



C<i,<I>k + C**k 



(16) 



with C$ = aa + (I — a)c and C$ = (1 — a)b + ac, where, 
in the limit w ^ 0, a describes the area of the cold 
region, a — $k, and 1 — a the area of the hot region, 
1 — a — ^jj. 4'k. AH the sums are normalized with respect 
to the number of fc-points of the BZ. In the limit w 
we get two different lifetimes in the cold and in the hot 
regions: 



lim Tk = 



— ^- in cold regions, , ^ 

— 7-4 — rr in hot regions. 



We consider the following temperature dependences of 
the scattering amplitudes in the scattering matrix: 



a{T) ^ aT^- b{T) = 6; c{T) ^ cT, 



(18) 



with the temperature independent parameters a, b and 
c having proper dimensions given by [?ia/A:^] = 1/eV, 
[hb\ = eV and [hc/ks] = 1. Indeed, in the cold region 
the scattering is weak and a Fermi liquid behavior with 
a temperature dependence of the scattering matrix 
Ck,k' (with both k, k' inside the cold region) is a reason- 
able assumption. On the other hand, in the hot region 
the scattering is strong and, as suggested by ARPES, 
the states are almost incoherent; we consider Ck,k' (with 
both k, k' inside the hot region) as temperature indepen- 
dent. Finally, we introduce a coupling between the hot 
and cold regions. The inter-patches elements of Ck.k' 
(with k inside the cold region and k' inside the hot re- 
gion and vice-versa) are considered to be temperature 
dependent with a linear in T dependence, and this is a 
key assumption in our model to obtain the linear temper- 
ature dependence of the resistivity, as shown in the next 
section. The important consequence of our assumptions 
for Ck,k' , and in particular the introduction of the inter- 
patches scattering with a linear temperature dependence, 
is that at low temperature the scattering amplitude is 
non- Fermi liquid at any k-point of the BZ. Indeed in the 
cold region we have 1/tc = aaT^ + (1 — a)cT, while in 
the hot region l/r/j = (1 ~ Q^)^ + o:cT. The Fermi liq- 
uid behavior can be recovered in the cold region when 



the area of the hot region tends to zero (i.e. a = 1). 
The non Fermi liquid behavior of the scattering ampli- 
tude in the cold region at low temperatureJs supported 
by the ARPES experiments of Valla et aZBj, as already 
discussed. 

To solve the BE, it is necessary to find how the opera- 
tor K"^ acts on an arbitrary velocity, as can be seen in 
Eq. (^. Because of the symmetry properties of the scat- 
tering matrix Ck,k' and of the quasiparticle velocity Atj-, 
when summing over the whole BZ, Ck.k'Vk' = (O. 
Using these properties, we find that acts on an ar- 
bitrary velocity just by inserting a scattering time Tk, 
that is momentum dependent, in front of the velocity, 
X^k' -^k k'''^k' — ■^k^'k- This allows us to solve exactly 
the linearized BE and to obtain formulas for the several 
transport properties we are interested in. The deviation 
(jk from the equilibrium distribution function is = 
ei?TkVk-n(— 9/^/9ek), where n is the direction of the ex- 
ternal electric field E = Eh. The momentum dependence 
of 5k is therefore given by gw ^ ■ n/(C$$k + C^'^k) 
and hence in the hot regions g^ is strongly suppressed 
while in the cold regions g^ has sizeable values. Hlubina 
and Rice, solving the BE by a variational approach in the 
case of hot spots generated by antiferromagnetic fluctu- 
ations coupled to fermions, obtained a g^ with similar 
properties, i.e. a depopulation of the hot regional^ (see 
Section IV for a discussion). Finally, as regard the sin- 
gle particle electronic properties, we recall that optimally 
and overdoped cuprates display very similar shape and 
momentum dependence of the Fermi surface and of the 
electronic band dispersion. Therefore we consider Bi2212 
as typical for the FS and band structure of cuprates. For 
this material several ARPES measurements are available. 
We consider a tight-binding model for the band structure 
of Bi2212 with hopping up to the fifth nearest neighbors. 
We are using the following fit for the enewf of the quasi- 
particles e(k) obtained by Norman et alE3. 



e(^kx, ky) — ^ ^ CiTjii^kx , ky 



(19) 



i=l 



with the values of the coefficients q and the func- 
tions rji given by ci = O.lSOSeF, C2 — — 0.5951ey, 
C3 = 0.1636ey, C4 = -O.OSlQeF, cg = -O.llUeV, 
C6 = 0.0510ey and f]i — I, r]2 ^ 0.5[cos{kx) + cos{ky)], 
ri3 — cos{kx) cos{ky) , rji — 0.5[cos(2A:3;) -I- cos(2fcy)], 
775 — 0.5[cos(2A:a;) cos(A:y) -I- cos(2fcy) cos(fc^)], 775 = 
cos(2fc2;) cos(2fcj,). These parameters are appropriate for 
the band structure of _Bz2212 compounds at optimum 
doping, thus giving an open FS and a van Hove singu- 
larity (VHS) slightly below the Fermi level. The value 
of — ci = — 0.1305ey is fixed to have the proper 

distance of the Fermi level from the VHS {ep — ^vhs = 
35meV as determined experimentally), and corresponds 
to the optimum doping 5 = 0.17. The bandwidth of the 
dispersion given in Eq. (p^ is lAeV. 
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III. TRANSPORT PROPERTIES 

A. Resistivity 

Theoretical results. The leading order contribution to 
the DC-conductivity is given by the first term in Eq. (4), 
thus A^\_, is replaced by K^]^, and inserted into Eq. (7). 

The operator K^^,, applied on a velocity Wk gives WkTk, 
as already discussed in Section |l[ Therefore, the DC- 
conductivity is given by: 



2e- 



E 



9ek 



(20) 



In the limit T — > we can obtain a finite contribution 
to fj^^ depending on the width w and on the angle 9 
(determining the area of the hot region). As Tk given 
in Eq. ( p^ ) doesn't diverge in the case w ^ (because 
^'k 7^ even inside the cold region), we get a finite con- 
tribution to a^^ and a residual resistivity at zero temper- 
ature. Note that the residual resistivity is determined by 
the width w and 6(1 — a), thus the T-independent part 
in the coefficient . In the limit w ^ we separate two 
regions in the BZ and we can consider a hot and a cold 
average velocity, and Vc respectively. In the limit of 
low T and w — > 0, the DC-conductivity is given by 



lim 



= [vlTcNciep) + vlThNhiep)] ■ (21) 



The quantity Nc^ep) is the density of states at the Fermi 
level in the cold region, iVc(e) = X]k'^(^k — e)^'k, while 
Nhiip) is the same quantity in the hot region, with 
Nhle) = J^k^i*^^ ~ ^)*k- Fig. I shows the density of 
states in the hot (upper panel) and in the cold region 
(lower panel) as a function of the size of the hot region 
(determined by the angle 9) in the limit i/j — s- 0, obtained 
using the band dispersion given in Eq. (p^. Eq. (21) 
shows that cold and hot regions are wired in parallel, 
while each region (cold/hot) is wired in series with the 
transition region (in the case w — 0), as shown by the 
structure of Tc and given in Eq. (p^). In the lim- 
its mentioned above the temperature dependence of the 
resistivity up to second order in T is given by: 



lim = 



1 (l-a)c 



T 



(22) 



1 



a c 



T 



As both lifetimes Tc and for w = diverge in the limit 
T ~^ 0, the resistivity has no zero temperature offset, 
p{T — > 0) = 0. However an offset can be obtained by a 
finite width w as discussed above. 
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FIG. 3. The overall shape of the hot density of states A^h(e) 
(upper panel), doesn't change with 6, while the magnitude at 
the Fermi level changes. In the case of the cold density of 
states Nc(eF) (lower panel), both the overall shape and the 
magnitude at the Fermi level changes. The position of the 
Fermi level Ep, corresponding to optimum doping, is also 
reported. Note the changement of sign in the derivative of 
Nciep) for e = 12.9°. 

In our two-patch model we can associate the dop- 
ing variation of the electronic properties of the cuprates 
mainly with the variation of the area of the hot region 
and hence with the angle 9. Indeed, ARPES experiments 
on Bi2212 show that the region of the BZ where the 
spectral function is broad and no quasiparticle peaks are 
detectable increases as the doping is reduced. A small 
segment of (quasiparticle) FS is observed approaching 
the metal-insulator transitionEl. Therefore, we can trans- 
late this behavior of the FS and of the spectral function 
saying that when the doping is decreased, the angle 9 
increases (i.e., the size of the hot region increases). Fig. 
^ shows the change in resistivity with increasing angle 
9; we obtain that the residual resistivity increases with 
increasing angle 9. 
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FIG. 4. Resistivity p(r) for different sizes of the hot region 
(i.e. different ff) as a function of temperature. The variation 
of the hot region changes the offset of the resistivity. It can 
be seen that increasing the area of the hot region increases 
the residual resistivity po and change slightly the power law 
behavior p(T) ~ T^. The following parameters are used: 
a = 48, 6 = 2, c = 7 and w = 0.3. 9 changes between 10° and 
30°. 

As can be seen in Eq. (^2|) the hot/cold scattering 
term c is the most important quantity in determining 
the slope of the resistivity respect to the other scattering 
amplitude a and h. The variation of the slope with c is 
reported in Fig. ^. Note that a variation in c doesn't 
change the residual resistivity. 
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FIG. 5. The amplitude of the scattering between the cold 
and the hot region, described by c, determines the slope of 
the resistivity. Note that po doesn't change. Increasing the 
inter-patch scattering determines an increasing in the slope 
of the DC-resistivity. The following parameters are used: 
a — 48, b = 2, 9 = 20° and w — 0.3. c changes between 
5 and 9. 

Moreover, the slope of the resistivity is controlled by 
the single particle properties of the cold region, in par- 
ticular by the cold density of states Nclep) and by the 



Fermi velocity v"^, and by the area of the cold region a, 
as obtained in Eq. (p^. 

On the other hand, the influence of the prefactor (1 — 
a)/ {v'^Nc{eF)) on the slope of the resistivity is small in 
comparison with the influence of c, as can be seen in Fig. 
^ and Fig. because of compensation effects: VcNc{eF) 
is roughly 9 independent and also the ratio (1 — ol)/vc is 
expected to vary smoothly with 9. 

Comparison with experiments. The linearity of the re- 
sistivity up to veryJiigh temperatures was found exper- 
imentally (see, e.g.l3). Experiments show a quasi-linear 
temperature dependence of the resistivity p^^{T) (x T'', 
where in Reffl it is shown that 7 increases with doping. 
For optimally doped systems 7 = 1, while for overdoped 
systems a 7 > 1 (7 ~ 1.5) is observed, supporting a 
gradual recovering of the Fermi liquid properties. The 
doping dependence of the exponent 7 can be compared 
with our results in Fig. |^, showing that small values of 
9 {9 < 15°) give a resistivity p{T) ~ with 7 > 1. 
Moreover the residual resistivity po — p{T = 0) increases 
with underdoping and this is in agreement with the trend 
obtained in Fig. |^. Therefore, the two-patch model 
provides the possibility to change the doping mainly by 
changing the angle 9, besides considering the doping de- 
pendence oi ep and WF(k). Note that the range of thg 
change in p with changing doping is observed in Ref.Q 
as well. Fig. 1 of Ref.Ll shows that for undcrdopcd 
Bi22Q\ (with a La concentration of a; = 0.66) the offset 
Po — 150/ir2cm, while for the overdoped compound (with 
X — 0.24) Po ~ 40pilcm. The increase in the slope and 
in the offset of the resistivity is also observed in Y12?> 
as the doping is reducedlJ: the offset changes between 
250pncm and 20pf7cm (Fig. 1(a) of Ref.El). The same 
range of variation is given in our results of Fig. ^ con- 
sidering the range 10° < 9 < 25°. Within the described 
change in 6* a big variation of doping can be described. 

In the following paragraph we report a simple explicit 
comparison between the transport properties here evalu- 
ated and the experimental data via a fitting procedure. 
The comparison here presented is limited to the Bi-hased 
cuprates {Bi2212 and Bi2201) for which several ARPES 
and transport experimental data are available. 

The strategy of fitting the experimental data is the fol- 
lowing. We introduce a temperature scale To upon which 
the term linear in T in the scattering rate in the cold re- 
gion dominates. This guarantees us that the resistivity 
is linear up to this temperature To. We choose a tem- 
perature To that has a value of « lOOOif . The relations 
a — (1 — q;)c/(q;To) and b — acTo/{l — a) allows us to ob- 
tain first values of a and b for given To and c. The width 
of the transition region w is fixed mainly by the mag- 
netoresistance as shown in section III-D. Starting with 
different angles 9 {i.e. different sizes of the hot region), 
we try to get a good fit of Hall-angle data as shown in 
section III-C. Then we try to fit the Hall-angle and the 
slope of the resistivity for different combination of 9 and 
c. It turns out that the angle 9 « 20° and the value 
c « 6.5 gives a good fit of the Hall-angle and of the slope 
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of the resistivity. After that, we change the scattering 
amphtude b in the hot region. Increasing b leads to an 
higher resistivity and thus the parameter b ahows us to 
adjust the ofFset of the resistivity. In this manner we can 
fix the parameters w,9,c,b. The last parameter that is 
to be fixed is d. The freedom for the parameter d is not 
big, in order to maintain the linearity of the resistivity. 
Thus the initial condition we have chosen for d can re- 
main valid. We obtain a reasonable fit for the different 
transport quantities with the following values of the pa- 
rameters: w = 0.20, a = 60, & = 2.1, c = 7.0 and 9 = 20°. 
In Fig.^ we report the resistivity as a function of temper- 
ature evaluated by Eq. (|2^) with the set of parameters 
given above and we compare our results with the resis- 
tivity measured in JBd2212 at optimum doping as given 
in Ref.El and in Ref.c3. 
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FIG. 6. Fit of the temperature dependence of the resistivity 
for the set of parameters given in the text. The experimental 
data are from an optimally doped Bi2212 samples from Ref. 
[1] and from Ref. [43]. 



B. Thermoelectric power 

Theoretical results. To evaluate the thermoelec- 
tricpower (TEP), which is the quantity measured in 
experiments, we first compute the longitudinal ther- 
mopower given in Eq. (^). The thermopower has 
almost the same expression as the DC-conductivity be- 
sides the extra factor of and an extra minus in the 
sum, and it is given by 



S*^^ = -2, 



(23) 



A Sommerfeld-expansion of Eq. (^ ) , in the limit of w 
and at low temperature, gives 



with the first derivatives of the density of states N'^^^^ep) 
evaluated at the Fermi level in the cold and in the hot 
regions. We found that S*^^ has a smooth dependence 
on w and hence the limit mentioned above is valid also 
for small values of the width w. As in the case of con- 
ductivity, because Tc » and Vc >> Vh, the main 
contribution to the thermopower comes from the cold re- 
gion. It is interesting that the Fermi level is below the 
peak in the cold density of states for 9 > 10° (see Fig. 
(lower panel)). The peak in iVc(e) is exactly at the 
Fermi level for 9 « 10°. The TEP, as measured in the 
experiments mentioned above, is defined as 

TEP=j;^, (25) 

where a^^ is given in Eq. (|20|). Again the cold region 
has the main influence. The expression of the TEP in 
the limit of w = and at low temperatures is given by 



TEP = kiT 



3e 



3e 



-kiT 



N4eF) + [{vlr,,)/{v^,r.)]N^,{ep) 
N'Aep) 



(26) 



where the last relation is obtained in the limit Tc >> r/j 
and Vc » Vh- Therefore the ratio between the derivative 
of the density of states and the density of states evaluated 
at the Fermi level in the cold region is the main quantity 
that determines the slope of the TEP at low tempera- 
ture. Increasing the temperature above 200i^r, we have 
verified that the contribution of the hot region to the 
TEP becomes comparable to the contribution from the 
cold one. 
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FIG. 7. Temperature dependence of the thermoelectric 
power for different angles 6. The following parameters are 
used: a = 48, fe = 2, c = 7 and w = 0.3. 6 changes between 
10° and 30°. 

Note that in the case of the TEP the offset doesn't 
change varying w. We found that the TEP is the most 
sensitive quantity to a variajtion of the electronic struc- 
ture, in agreement with Rcf.H. This is done by changing 
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the hopping parameters Ci and C2 in Eq. (p^). Fig. |^ 
shows the TEP obtained from the two-patch model for 
different angles 9. For 6* > 10° the TEP is negative and 
the slope of the TEP at low temperature changes slightly 
with different angles 9 > 10°. We found that, as pre- 
dicted in the limit given in Eq. (26), the TEP is almost 
not affected by other parameters than 9. Nonetheless we 
got the largest effect on it from the scattering in the hot 
region. Again the slope doesn't change so much and the 
offset changes also slightly when the hot-scattering b is 
changed of a factor 5, as shown in Fig. M. 
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FIG. 8. Scattering in the hot region slightly affect the slope 
of the TEP and it has a weak influence on its offset. An in- 
crease in the hot-region scattering lowers the TEP and makes 
it even more negative. The following parameters are used: 
a = 48, c = 7, 9 — 20° and w — 0.3. b changes between 0.5 
and 2.5. 
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FIG. 9. Comparison between the TEP evaluated by the 
two-patch model and the thermoelectric data given by Mcin- 
tosh et al. in Ref. [3] (material: Bi2201 at optimum doping 
(x=0.44) and overdoping (x=0.2)) and by Obertelli et al. in 
Ref. [46]. (material: Bi2212 at optimum doping {S = 0.0) 
and overdoping (5 = 0.077)). 

We attribute the quantitative discrepancy between the 
TEP measured in optimally and underdoped cuprates to 
the fact that the TEP reflects the properties of the exci- 
tations away from the FS and hence the frequency depen- 
dence of the scattering aipelitude becomes important, as 
shown by DMFT studiesc2l. We deserve the inclusion of 
the frequency dependence of the scattering amplitudes 
for future investigation. 

C. Hall angle 



Comparison with experiments. In Fig.^ we report the 
TEP evaluated by Eq. ( p3| ) with the same set of pa- 
rameters used for the resistivity {w = 0.20, a — 60, 
b = 2.1, c = 7.0, 9 = 20°) and we compare our results 
with the TEP measured in Bi220l at optimum and over 
doping as given in Ref.B and with the TEP measuredpLn 
Bi2212 at optimum and overdoping as given in Ref.c3. 
The TEP obtained with our two-patch model is only in 
qualitative agreement with the TEP observed by exper- 
iments. In particular only the overdoped Bi2201 has a 
TEP with a temperature dependence close to the one 
evaluated within the two-patch model. Note also that 
the TEP measured in overdoped T/2201 cuprates has al- 
most the same temperature dependence and magnitude 
of overdoped Sz2201 (and hence close to thejXEP of the 
two-patch model) , as shown in Fig. 1 of Ref.Ej. 



Theoretical results. The leading order contribution to 
the Hall-conductivity cr^*' is given by the second term in 
the expansion of the operator A^^ respect to the (weak) 
magnetic field, A^^ — —K^^MbK"^ . Inserting this 
term into Eq. (^, we get the Hall-conductivity a^^ . The 
bending term in A^^ is 

(vk X B)Vk = B{vldk^ - v^dkj, (27) 

that arises from the operator Mb considering a magnetic 
field perpendicular to the CuOi planes. In this case, the 
formula for the Hall-conductivity is the following 

k 

(28) 

The partial 
derivatives of the relaxation time Tk enter in Eq. ( |2^ ) 
and are given by dTi^/dk^ = —r^id^k) / (dkx) [C$ — C*] 
and dri^/dky — — T^(9<i>k)/(9fcjy) [C<i. — C^] respectively 
(note that d^-k/dkx — —d'^^^/dkx). The final expression 
for the Hall-conductivity is 
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k 



;^^k 



(29) 



+ 
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he 



k 



This expression contains two different powers of the scat- 
tering time, one oc and another oc r^. In the hniit 
w — > we can get some further insight in the problem of 
the HaU-conductivity. In the low temperature limit the 
sum over k is restricted on the FS. The HaU-conductivity 
contains in this limit derivatives of step functions and 
hence (5-functions are generated. As a consequence we 
get 8 points on the FS that contribute to the second 
and third term in Eq. (29). The temperature depen- 
dence of the relaxation time tq evaluated in the tran- 
sition region hot/cold on the FS can be written, using 
in this region <i>k = ^'k = 1/2 and using Eq. (|lq), as 
To = 1/(2C<E, + 2C*) = l/(co + ciT + C2T^). In this limit 
Eq. (Efl) can be written in a more compact form as 



{Vy{A)VyiB) [Vy{B) 



Vy{A)]}[q 
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where a^^ is the first term in Eq. (p9|), A and B label the 
contributing points on the FS in the first quadrant of the 
BZ and q is either tan 9 or cot 6 depending on 9 respect to 
the value 6 ~ 10°. Note that we get an offset to a^^ from 
the last terms even in the limit w ~^ 0. Increasing w, also 



the first term in (29) contributes to the offset due to the 
same reasons given for the resistivity. The first term in 
( p9| ) contains the customer contribution proportional to 
T^while the second term contains higher powers in rt and 
can originate deviations from the behavior depending 
from the temperature and the scattering amplitudes. 

The cotangent of the Hall-angle is defined as the ratio 
between the direct conductivity and the Hall conductiv- 
ity, 



cot 9 HiT) = —. 



(30) 



A plot of cot 9h vs. (see FigJTo|) shows that the 
cotangent of the Hall-angle has a temperature depen- 
dence oc in the high temperature regime (T > 200K) , 
where we observe almost a straight line. 
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FIG. 10. The size of the hot region determines the offset 
of the Hall-angle. In the case of = 30°, which is equivalent 
to 3/4 of the first BZ with hot character, the offset in the 
Hall angle is close to 50. Note that we plot cotOn vs. T^. 
The following parameters are used a = 48,fe = 2,c = 7 and 
w = 0.3. 6 changes between 10° and 30°. The magnetic field 
is B = IT. 

It can also be seen in Fig. ^ that the range of temper- 
ature where cot 9h ~ increases with decreasing 9 (i.e. 
increasing doping). Konstantinovic et al.a pointed out 
the strong influence of the anisotropy of the FS on the 
Hall-angle. We reproduced this observation by changing 
the hopping parameters ci and C2 in Eq. ( p^ ) but we 
found that the effect of this change has even bigger ef- 
fects on the TEP, as already discussed. Among the other 
parameters we found the strongest influence on cot 9h (T) 
from the inter-patch scattering amplitude c. The effect 
of c on the Hall-angle is shown in Fig. pd] . 
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FIG. 11. Effect of the hot/cold inter-patch scattering on 
the Hall-angle. It turns out that c has an effect on the temper- 
ature dependence of the Hall-angle, but no significant effect 
on the offset of the Hall-angle contrary to 6. The following 
parameters are used: d — 48, b — 2, = 20° and w = 0.3. c 
changes between 5 and 9. The magnetic field is B = IT. 
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It can be seen that this parameter changes the tem- 
perature range where cot Oh ~ only slightly and the 
offset is independent on c. The larger the value of c is 
the smaller the deviation of the Hall-angle from a straight 
line is, which means that increasing the cold/hot scatter- 
ing has opposite effect as decreasing the angle 9. On the 
other hand changing 6 (i.e. doping) effects the offset of 
the Hall-angle. In the framework of the two-patch model 
we obtain a decreasing offset of the cotangent of the Hall- 
angle with increasing doping and an increasing value of 
the power law coefficient 7 (cotOfj oc T^) with increasing 
doping. 

Comparison with experiments. As already discussed, 
the range of temperature where cot6'_ff ~ increases 
with decreasing 6 (i.e. increasing doping V|(see Sig. 10) 



which is contrary to the results |given inO andEl, but in 
agreement with the results of Refo. Indeed a careful ex- 
amination of various set of experimental data for cot 6h 
vs. indicates that the exact dependence does not 
extend over the whole temperature range, and jieviations 
are observed for T < 300K. As shown in Ref.EJ, e.g., in 
y 123 as Tc is lowered by underdoping, the range of the 
dependence of cot Oh moves toward higher temperatures 
(for Tc = 90K, the deviation start at T = lOOK, while 
for Tc = 40iC, the deviation start at T = 170K). On the 
other hand changing 6 (i.e. doping) effects the offset of 
the HaU-angle which was found inQ. Indeed in Fig. 3(b) 
of Ref.Q the change in the offset of the cotangent of the 
Hall angle between underdoped (x=0.66) and overdoped 
(x==0.24) Bz2201 compounds measured at T ~ 280K 
(T^ ~ 8 • lO^K^) is of the order of 1300, a value which is 
compatible with our finding of Fig. |l^ (~ 2000) consid- 
ering again 10° < 9 < 25° as in the case of the resistivity 
discussed in Section III-A. In the case of yi23 the slope 
of cot9H vs. increases when the doping is reduced 
(almost 30% of variation, when Tc is reduced from 90A' 
for optimally doped to 40Ar for an underdoped yi23, 
while in this case the offset increases-only slightly (see 
Fig. 1(e), Fig. 3 and Fig. 4 of Ref.El). Note that for 
yi23 a non monotonic behavior of the slope of cot 6h 
vs. as a function of doping is obseriied in a differ- 
ent set of measurements reported in Ref.EJ. This can be 
understood remembering that we increase the hot region 
increasing 9 (i.e. decreasing doping). In the framework 
of the two-patch model we obtain a decreasing offset of 
the cotangent of the Hall-angle with increasing doping 
and an increasing value of the power law coefficient 7 
{cot 9h oc T'^) with increasing doping. 

In Fig.|l^ we report the cotangent of the Hall angle 
evaluated by Eq. (|3^) with the same set of parameters 
used_for the resistivity and the TEP {w = 0.20, a 
60, b = 2.1, c = 7.0, 9 = 20°) and we compare our 
results with the cotangent of the Hall angle measured 
in _Bz2212 at optimum doping as given in Ref.El |and in 
_Bi2201 again at optimum doping as given in Ref.lZl. The 
agreement with the data obtained in Bi2212 is good and 
also a qualitative agreement is obtained if we compare 
the results obtained with the two-patch model and the 



band structure of _Bi2212 with the experimental data for 
Bi2201 at optimum doping. 
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FIG. 12. Comparison between the experimental data for 
the Hall angle and the Hall angle obtained from our two-patch 
model using the set of parameters given in the text. The 
experimental data that are shown in the plot correspond to 
5*2212 at optimum doping from Ref. [6] and to _Bi2201 at 
optimum doping from Ref. [7]. 



D. Magnetoresistance 

Theoretical results. The magnetoresistance (MR) is 
defined as the ratio of the variation in the resistivity 
in presence of a magnetic field to the resistivity with- 
out magnetic field. In the transverse geometry, with the 
magnetic field applied perpendicular to the Cu02 planes 
and the current measured parallel to the planes, the MR 
is given by 



MR- 



Ap^^(B) 
p^^(O) 



Acr^^(B) 
cr^^(O) 



tan^ 9h. 



(31) 



In this geometry, the first extra contribution to the 
DC-conductivity, Aa^^{B) is achieved by the third 
term in Eq. (jj), thus we have to insert — 
K-^MbK^^MbK^^ into Eq. (0). We get an expres- 
sion for the correction to the conductivity given by 
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— In this case it is necessary to consider partial deriva- 
tives of with respect to or ky, which are 
given by dr^/dk^ = -2T^9$k/9/ca; [(7$ — C*] and 
drl/dk^ = -3r4(a$k)/(afc^) [CcE, - C^]. Finally we ob- 
tain Aa^'^'iB) as: 
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We get three contributions to Aa^^(i?) with different 
powers of Tk, c!c T^, oc and oc t^. This formula is in- 
serted in Eq. (|3|) together with the Hall-angle computed 
in the previous subsection. We are now able to compute 
the MR for different parameters. As shown in Fig. a 
change in the angle 9 (and hence in doping) has a size- 
able effect on the MR and decreasing the angle 9 (i.e. 
increasing doping) increases the MR in our model. 
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The effect of the transition width w is studied in Fig. 
It can be seen that this quantity becomes more and 
more important the smaller it becomes. A big difference 
in the MR can be observed between w — 0.1 and w — 
0.2, which is in agreement with the results for the MR 
evaluated within^ cold spots model by Zheleznyak et al. 
reported in Refa. 

Comparison with experiments. Ando et al. reported 
that the orbital MR increases with increasing doping, as 
shown in Fig. 5 of RefE. (The orbital contribution to 
the MR can be obtained by the transverse component of 
MR subtracting the longitudinal component, eliminating 
in this way the contribution of the spins to the MR.) This 
experimental result is derived in our model. In F ig, we 
report the magnetoresistance evaluated by Eq. ( |3l| ) with 
the same set of parameters used for the other transport 
properties above discussed {w — 0.20, d — 60, b — 2.1, 
c = 7.0, 9 = 20°) together with B = IT and we compare 
our results with the orbital magnetoresistance measured 
in Bi2201 at ogtimum doping and slightly underdoping 
as given in Ref.Ll. (MR data for Bi2212 are not yet avail- 
able to our knowledge). For the high doping level and 
range of temperature here considered the longitudinal 
MR is an order of magnitude smaller than the transverse 
MR and its contribution to the orbital MR is therefore 
small. 
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FIG. 13. The effect of variation in the angle 6 on the MR. 
The MR decreases increasing the area of the hot region. In 
the limit where the whole BZ is hot {6 = 45°) the MR is zero. 
The parameters are a = 48, 6 = 2, c = 7 and w — 0.3. 9 
changes between 10° and 30°. The magnetic field is _B = IT. 
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FIG. 14. The MR diverges in the limit w ^ 0. The influ- 
ence of w on the MR is very strong, so the MR allows us to fix 
the parameter w quite well as other quantities don't depend 
strongly on w. We used the parameters a = 48, 6 = 2, c = 7 
and 8 — 20° in the plot, w varies between 0.1 and 0.5. The 
magnetic field is i? = IT. 
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FIG. 15. Comparison between the magnetoresistance ob- 
tained from the two-patch model and the orbital magnetore- 
sistance of an optimally doped (a; = 0.44) and sHghtly under- 
doped {x = 0.57) Bj2201 from Ando et at [7]. 

While the order of magnitude and the qualitative tem- 
perature dependence of the MR evaluated with the two- 
patch model agree with the MR data for 5^2201, a quan- 
titative discrepancy is observed taking fixed the set of pa- 
rameters we used to evaluate the previous transport prop- 
erties. On the other hand a change in w from w — 0.2 
to w = 0.1 can increase the MR curve shown in Fig. |l^, 
leading to a better agreement with the low temperature 
MR data for Bi2201. Note that the width of the transi- 
tion region w is associated to the properties of the elec- 
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tronic scattering and its value can be material dependent, 
even in presence of a similar FS and band structure. A 
temperature dependence of the coefficient w seems neces- 
sary to improve the fit of the MR data. MR data for op- 
timally doped i?z2212 could permit a more careful com- 
parison with the prevision of the two-patch model here 
presented. 



E. Thermal Hall conductivity 

Theoretical results. The thermal-Hall conductivity is 
evaluated using Eq. (||). In this case the operator A'^ 



has the form K ^MbK ^. Replacing 
(p9|) gives us the result for k^^. 



in Eq. 
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suppress in a sizeable way. The same behavior is 
observed considering the role of the hot patches, changing 
the scattering amplitude b. 
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In Fig. |l^ the thermal-Hall conductivity is reported as 
a function of the temperature. 



FIG. 16. The thermal-Hall conductivity (like electrical 
transport properties) is mostly influenced by the cold region. 
We use the following parameters: d — 4S, b — 2, c — 7 and 
w = 0.3. 9 changes between 10° and 30°. 



FIG. 17. A change in c shifts the thermal Hall conduc- 
tivity: an increasing of the inter-patch scattering leads to a 
decreasing of the thermal-Hall conductivity. The following 
(34^arameters are used: d = 48, b = 2, w — 0.3 and 6 — 20°. c 
changes between 5 and 9. 

Considering the qualitative correspondence between 
decreasing doping and increasing 9 discussed above, the 
two-patch model suggests that in underdoped cuprates 
the thermal-Hall conductivity should be strongly sup- 
pressed with respect to the one measured in optimally 
and overdoped cuprates. This conclusion is also sup- 
ported by the fact that the electronic scattering in the 
hot region (and hence b) should increase in the under- 
doped regime due to the proximity to the antiferromag- 
netic phase. 

Comparison with experiments. In Fig. ^ we report 
the thermal-Hall conductivity evaluated by Eq. (34 ) with 
the same set of parameters used for the other transport 
properties above discussed {w = 0.20, a = 60, 6 = 2.1, 
c = 7.0, 9 = 20°). Again, experimental data for k^^ 
for the _Bi-based cuprates are to our knowledge not yet 
available. The temperature dependence of k^^ agrees 
quite well as regard the magnitude with the data for k^^ 
measured for an optimally doped YBCOeI, even if this 
material has some differences in the band structure and 
FS respect to the Bi-hased cuprates. 



The various curves correspond to different values of the 
angle 9 and we obtain a large increase of k^^ as the angle 
9 is decreased. As in the case of the other transport 
properties, the cold patches has the main influence in 
determining k^"^, even if the contribution to k^^ from 
the hot patches is sizeable as in the case of the TEP. 
The role of the inter-patches coupling in is studied 
changing the scattering amplitude c and the results are 
reported in Fig. |l^, showing that increasing c tends to 
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FIG. 18. Thermal-Hall conductivity k^^ vs. temperature 
obtained from the two-patch model with the same set of pa- 
rameters given in the text. The temperature dependence of 
hi^^ agrees quite well as regard the magnitude with the data 
for measured for an optimally doped YBCO [8]. 



IV. CONCLUSION AND DISCUSSION 

Normal state transport properties of cuprate supercon- 
ductors have been studied using a semi-classical approach 
based on the hnearized Boltzmann equation (BE). The 
probabihty of scattering between the electrons of the con- 
duction band and an effective collective mode is assigned 
by a scattering matrix defined on a Brillouin zone (BZ) 
divided in two kind of patches. One kind of patches, 
centered in the M points of the BZ, contains hot states, 
which are strongly interacting, and is characterized by a 
low Fermi velocity and an high local density of states. 
The other kind of patches, centered in the nodal points 
along the BZ diagonals, contains cold states, which are 
weakly interacting, and is characterized by an high Fermi 
velocity and a low local density of states. In the multi- 
patch model the scattering matrix is assumed to be a 
sum of separable terms (in k and k') having coefficients 
with a temperature dependence which can be non con- 
ventional, according to the properties of the scattering. 
In the case of the two-patch model, the scattering ma- 
trix is a sum of a term describing the scattering between 
the electrons inside the hot patches with a large ampli- 
tude and temperature independent, a term describing the 
scattering inside the cold patches with a smaller ampli- 
tude and dependent on the temperature as T^, and a 
term including in a symmetric way the inter-patch scat- 
tering, with a linear temperature dependence. With this 
phenomenological choice of the scattering matrix, the BE 
is exactly solvable, and all the transport properties (at 
least in the weak field regime) can be evaluated. The 
resulting scattering amplitude 1/rk is strongly momen- 
tum dependent (in particular along the Fermi surface) 
and the low temperature behavior is always non Fermi 
liquid, with a linear temperature dependence in the cold 



patches and a constant in the hot patches, as suggested 
by recent ARPES experiments performed on Bi221'£3. 
The deviation of the distribution function from the equi- 
librium is strongly suppressed in the hot patch because 
of the strong scattering, while in the cold patch it has a 
sizeable value, because of the weaker scattering, giving 
the largest contribution to transport. 

The two-patch model here introduced has similarities 
as well differences with the model of Hlubina and Ricell^. 
Hlubina and Rice consider a model where the fermions 
are scattered by an antiferromagnetic spin fluctuation. 
The propagator of the spin fluctuations is peaked at the 
antiferromagnetic wave- vector Q = (tt; tt) and hence cou- 
ples mainly states around the M points, giving rise to 
hot spots, while the states around the nodal points are 
weakly coupled, giving rise to cold regions. The scat- 
tering matrix correspondent to this interaction has two 
different temperature regimes: (i) the low temperature 
regime, where the scattering between the hot states has 
a "v/T behavior, the scattering between the cold states 
has a Fermi liquid behavior and the scattering be- 
tween hot/cold states has again a behavior; {ii) the 
high temperature regime is instead similar to our two- 
patch model, having a constant scattering between the 
hot states, a quadratic (~ T'^) scattering between the 
cold states and a linear in T scattering between hot / cold 
states. Therefore, the model of Hlubina and Rice gives at 
low temperature a (Fermi liquid) temperature depen- 
dence of the resistivity, while our two-patch model gives a 
linear in T (non-Fermi liquid) behavior of the resistivity 
even at low temperature. 

An additive two-lifetime model, with similarities to our 
two-patch model, as been previously proposedcj. The re- 
sults obtained within this model are based on a peculiar 
Fermi surface, characterized by large flat regions paral- 
lel to the TM directions having a short relaxation time 
(hot regions) and small sharp corners around the nodal 
points along the TY{X) directions having a long relax- 
ation time (cold regions). Moreover the band structure 
considered in Ref.cll is such that the Fermi velocity is 
large in the hot region and small in the cold region. The 
authors suggest that this peculiar single particle proper- 
ties are typical for YBCO at optimum doping. ARPES 
experiments for Bi2212 and Bi2201 do not support this 
picture, and in particular the ratio between the Fermi ve- 
locity is the opposite, having large Fermi velocity in the 
cold region and small (almost undefined) Fermi velocity 
in the hot region. The flat region of Fermi surface ob- 
served in Bi-based comppunds is also much smaller than 
the one proposed in Rcf.o. This flat shape of the Fermi 
surface seems more appropriate for La-based underdoped 
cuprates, where stripe correlatipas can select preferred 
directions on the Fermi surfac^. The regions of the 
Fermi surface controlling the transport properties in the 
additive two-lifetime model and in our two-patch model 
are different. In particular the conductivity is controlled 
in the first case by the flat (hot) regions, while in our case 
is controlled by the curved (cold) regions, mainly because 
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of the completely different ratio between the Fermi veloc- 
ities. In the two-lifetime model the scattering amplitude 
in the hot region has a linear temperature dependence 
and a quadratic (Fermi liquid like) temperature depen- 
dence in the cold region. Therefore in this model the 
resistivity is linear in temperature, while the cotangent 
of the Hall angle is roughly quadratic in temperature, be- 
ing the Hall conductivity mainly controlled by the regions 
of Fermi surface with sizeable curvature (as the corners). 
In our approach both the resistivity and the Hall con- 
ductivity are mainly controlled by the cold regions of the 
Fermi surface, where, as already discussed, the scatter- 
ing amplitude has a non Fermi liquid character at low 
temperature. 

A further comparison Gari be done with the model pro- 
posed by loffe and Milliald, where the scattering ampli- 
tude is a sum of two terms, one temperature independent 
with an angle dependence along the Fermi surface with 
a deep (quadratic) minimum along the diagonal direc- 
tions {TY{X)) and the other Fermi liquid-like, with a 
quadratic temperature dependence. At low temperature 
this model gives a large scattering amplitude constant in 
temperature in the (large) hot region and a scattering 
amplitude with a Fermi liquid temperature depen- 
dence in the (small) cold region (cold spots). The linear 
resistivity is obtained in this model because the conduc- 
tivity is dominated by a small (cold) region which has a 
length proportional to the temperature, where the scat- 
tering is Fermi liquid-like. In our two-patch model the 
area (and hence the length) of the cold region is consid- 
ered to be temperature independent and a direct com- 
parison with the model of loffe and Millis is not possible. 
The form of the scattering amplitude proposed in this 
cold spjOts model is similar to the form proposed by Valla 
et alEB to describe the temperature and momentum de- 
pendence of the width of the (quasiparticle) peak in the 
ARPES spectra. On the other hand the experiments are 
consistent with a linear temperature dependence of the 
width around the zone diagonals. Interestingly, a quasi- 
particle peak is present in the ARPES spectra for all the 
Fermi wave- vector along the curved area of the Fermi sur- 
face, which is clearly not only a sharp corner, but a size- 
able fraction (almost the half) of the whole Fermi surface. 
Of course a direct comparison between the scattering am- 
plitude or the scattering matrix elements proposed in the 
various models and the ARPES line-shape is not possible 
and only some qualitative understanding can be obtained 
from ARPES experiments without a microscopic theory 
which is able to connect two-particle and single-particle 
properties. 

Another model for magnetotransport in cuprates has 
been recently proposed by Varma and Abrahamal3. This 
model combines the marginal Fermi liquid hypothesis for 
the inelastic scattering rate (linear in T), with the hy- 
pothesis that small-angle forward scattering is acting in 
the cuprates due to the scattering of the electrons of the 
Cu02 planes with the out of plane impurities; the for- 
ward scattering term in the scattering rate results tem- 



perature independent and strongly anisotropic. The total 
form of the scattering rate is again similar to the form 
proposed by Valla et a/.c2l to fit the ARPES spectra of 
Bi2212 as already discussed. Solving the BE, the au- 
thors show that forward scattering is responsible for a 
correction term in tanOn respect to the customary con- 
tribution and this term has the temperature dependence 
of the resistivity squared. In our approach the correction 
term corresponds to the second and third terms in Eq. 
(p9|), proportional to t^. On the other hand in our two- 
patch model we obtain a finite contribution from the first 
term, proportional to (the customary term), which is 
of the order or larger of the other terms. A detailed com- 
parison between the two approaches requires to consider 
a five-patch model as discussed in Section II. 

Our systematic analysis of normal state transport 
properties of cuprate superconductors, including resis- 
tivity, thermoelectric power. Hall conductivity, magne- 
toresistance and thermal-Hall conductivity, permits to 
understand which is the role of the patch geometry, of 
the Fermi surface and of the scattering matrix elements 
in determining the magnitude and the temperature de- 
pendence of the various electrical and thermal transport 
properties. In particular, the linear temperature depen- 
dence of the resistivity is associated to the inter-patch 
scattering, and its slope is determined by the amplitude 
of the inter-patch scattering but also by the single par- 
ticle properties in the cold patch, as the cold density of 
states evaluated at the Fermi level and the cold Fermi 
velocity. The Hall conductivity is also governed by the 
cold patch and the interplay between the various power 
law coming from the scattering time and from its par- 
tial derivatives gives the possibility to obtain a cotan- 
gent of the Hall angle with a T'^ behavior (with 7 > 1.5) 
in a range of temperature where the resistivity is lin- 
ear. The different power law behavior of resistivity and 
cotangent of the Hall angle indicates that the momentum 
dependence of the scattering time along the Fermi sur- 
face plays an important role and can originate a different 
characteristic scattering time for longitudinal and trans- 
verse transport. The magnetoresistance is mainly deter- 
mined at low temperature by the inter-patch region, and 
in particular its expression contains the second derivative 
of the scattering time. Therefore, the transition between 
the hot and the cold patches has to be smooth, in or- 
der to avoid a spurious divergence. Thermal transport 
is also considered in our systematic analysis. Thermo- 
electric power (TEP), as thermal-Hall conductivity, are 
mainly determined by the cold patch in the low tempera- 
ture regime, and it is interesting to note that the slope of 
the TEP is only controlled by the cold density of states 
and its derivative, giving a strong connection between the 
experimental slope of the TEP and the patch geometry. 
Increasing the temperature (T > 200K), the hot patch 
starts to contribute to the thermal properties. 

Finally, we present a tentative application of our two- 
patch model to the electrical and thermal transport prop- 
erties of optimally and overdoped iJz-based cuprates. We 
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use the electronic band structure and Fermi surface ob- 
tained by ARPES experiments and we obtain other in- 
formation by the temperature dependence of the ARPES 
Hneshape. The amphtudes of the intra-patch and inter- 
patch terms in the scattering matrix are fixed using the 
measurements of the resistivity, while the patch geometry 
is fixed by the TEP. Once all the parameters have been 
fixed, Hall conductivity, magnetoresistance and thermal- 
Hall conductivity are evaluated without any other as- 
sumptions and a reasonable agreement between our two- 
patch model and experimental values is found. In con- 
clusion, the two-patch model for the scattering process 
emerges as a minimal division of the Brillouin zone to 
account for the strong anisotropy of the eff'ective electron- 
electron interaction present in the cuprates, which is able 
to describe the several anomalous temperature depen- 
dences observed in the normal state transport proper- 
ties of cuprate superconductors. The description of the 
transport properties here presented can be improved in- 
creasing the number of patches, e.g. to N= 5 to include 
forward scattering processes. 

The functional form of the scattering operator can be 
realized in C-DMFT calculations. Work to see if a micro- 
scopic model such as a Hubbard model for some choice 
of parameters produces a temperature dependence close 
to our optimal fit of the data is currently under investi- 
gation. 



ACKNOWLEDGMENTS 

We appreciate valuable discussions with C. Castel- 
lani, M. Cieplak, M. Civelli, D. Drew and R. Raimondi. 
A. Perali acknowledges partial support from Fondazione 
" Angelo della Riccia" . M. Sindel is grateful to DAAD for 
partial support. 



E-mail addresses: perali@physics.rutgers.edu, mlsin- 

del@physi cs.rutRers.edu or k otIiarQphysics.rutRers.edu. 

Website: littp://www. htcs.org 
^ L. Forro, D. Mandrus, C. Kendziora, L. Mihaly, R. Reeder, 

Phys. Rev. B 42, 8704 (1990). 
^ K. Takagi, B. Batlogg, H.L. Kao, J. Kwo, R.J. Cava, J.J. 

Krajewski, W.F. Peck, Phys. Rev. Lett. 69, 2975 (1992). 
^G.C. Mcintosh, A.B. Kaiser, Phys. Rev. B 54, 12569 

(1996). 

" T.R. Chien, Z.Z. Wang, and N.P. Ong, Phys. Rev. Lett. 
67, 2088 (1991). 

5 A. , : , ~; 

Li, H. Raffy, Phys. Rev. B 62, 



Malinowski et ai, zo 
® Z. Konstantinovic, Z.Z 
11989 (2000). 

Y. Ando,T. Murayama, Phys. Rev. B 60, 6991 (1999) 



* Y. Zhang, N.P. Ong, Z.A. Xu, K. Krishana, R. Gagnon, L. 

Taillefer, Phys. Rev. Lett. 84, 2219 (2000). 
^ M.R. Norman et al, Nature 392, 1587 (1998). 
^° P.W. Anderson, Phys. Rev. Lett. 67, 2092 (1991). P.W. 

Anderson, The Theory of Superconductivity in the High- 

Tc Cuprates (Princeton University Press, Princeton, 1997), 

and references therein. 
" L.B. loffe, A.J. Millis, Phys. Rev. B 58, 11631 (1998). 

A.T. Zheleznyak, V.M. Yakovenko, H.D. Drew, I.I. Mazin, 

Phys. Rev. B 57, 3089 (1998). 
" R. Hlubina and T.M. Rice, Phys. Rev. B 51, 9253 (1995). 
^"^ A. A. Abrikosov, Fundamentals of the Theory of Metals 

(North-HoUand, Amsterdam, 1988). 

The hot /cold division of the BZ has also been used in the 
two-gap model to describe the pseudogap opening in under- 
doped cuprates induced by pair fluctuations as discussed in 
Ref.[a. 

A. Perah, C. Castellani, C. Di Castro, M. GriUi, E. Piegari 
and A.A. Varlamov, Phys. Rev. B 62, R9295 (2000). 
" H. Ding et ai, Phys. Rev. Lett. 76, 1533 (1996). 
J. Mesot et ai, Phys. Rev. B 63, 224516 (2001). 
J.C. Campuzano et ai, Phys. Rev. B 53, R14737 (1996). 
T. Valla et ai, Phys. Rev. Lett. 85, 828 (2000). 
D.S. MarshaU et al, Phys. Rev. Lett. 76, 4841 (1996). 
A. Carrington, A. P. Mackenzie, C.T. Lin, and J.R. Cooper, 
Phys. Rev. Lett. 69, 2855 (1992). 

C. Castellani, C. Di Castro, M. GriUi, Phys. Rev. Lett. 75, 
4650 (1995). 

A. Georges, G. Kotliar, W. Krauth and M.J. Rozenberg, 
Rev. Mod. Phys. 68, n. 1 (1996) and reference therein. 
^'^A. Schiller and K. Ingersent, Phys. Rev. Lett. 75, 113 
(1995). 

^'^ M. H. Hettler et el., Phys. Rev. B 58, 7475 (1998). 

Th. Maier, M. JarreU, Th. Pruschke, J. Keller, Phys. Rev. 
Lett. 85, 1524 (2000). 

C. Huscroft et el, Phys. Rev. Lett. 86, 139 (2001). 

A.I. Lichtenstein and M.I. Katsnelson, Phys. Rev. B 62, 
R9283 (2000). 

^° G. KotUar, S.Y. Savrasov, G. Palsson and G. Biroh, Phys. 
Rev. Lett. 87, 186401 (2001). 

G. M. Ehashberg, Soviet Phys. JETP 14, 886 (1962). 

D. Pines, Z. Phys. B 103, 129 (1997). 

F. Becca, M. Tarquini, M. GriUi, and C. Di Castro, Phys. 
Rev. B 54, 12443 (1996). 

A. Perah, C. Castellani, C. Di Castro, M. Grilh, Phys. Rev. 
B 54, 16216 (1996). 

G. Biroh and G. Kotha r, |cond-mat/0107108 



S. Onoda, M. Imada, :ond-mat/0108416, proceeding of 



the conference Spectroscopies of Novel Superconductors, 
Chicago (2001). 
■^^ D. Pines and P. Nozieres, The Theory of Quantum liquids. 
Vol. I, (W.A. Benjamin, Inc., New York, 1966). 
M.R. Norman, M. Randeria, H. Ding, J.C. Campuzano, 
Phys. Rev. B 52, 615 (1995). 

The zero average of Ck,k'v'k over the BZ is no longer valid 
when a generic many-patch model is considered. For in- 
stance, when forward scattering is included in Ck,k' , the 
scattering between different cold regions is not equivalent 
any more with the scattering inside a given cold region, 
and the average given above is in general different from 



17 



zero. The general solution for the multi-patch model given 
in Eqs. pj|]l^ , |l^ has to be considered. 

G. Palsson and G. Kothar, Phys. Rev. Lett. 80, 4775 
(1998). 

P. Xiong, G. Xiao, X. D. Wu, Phys. Rev. B 47, 5516 (1993). 

B. Wuyts et al, Phys. Rev. B 47, 5512 (1993). 

M.A. Quijada, D.B. Tanner, R.J. Kelley, M. Onellion, 
Physica C 235-240, 1123 (1994). 

A.T. Zheleznyak, V.M. Yakovenko, H.D. Drew, Phys. Rev. 

B 59, 207 (1999) and reference therein. 

A. Fujimori, A. Ino, T. Yoshida, T. Mizokawa, Z.-X. Shen, 

C. Kim, T. Kakeshita, H. Eisaki, S. Uchida, Proceedings 
of NATO Advanced Research Workshop Bled 2000 "Open 
Problems in Srongly Correlated Eletron Systems", edited 
by P. Prelovsek and A. Ramsak (Kluwer Academic). 
S.D. OberteUi, J.R. Cooper, J.L. Tallon, Phys. Rev. B 46, 
14928 (1992). 

C. Varma, E. Abrahams, Phys. Rev. Lett. 86, 4652 (2001). 



18 



